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Abstract: The standard argument why gravity is not renormalisable relies on direct pow- 
ercounting of Feynman graphs to estimate the degree of UV divergence. This analysis has 
in several (highly) supersymmetric examples be shown to overestimate divergences consid- 
erably. In these examples the main improvements arise from a conjectured duality between 
color and kinematics. In this paper we initiate the systematic study of quite general pow- 
ercounting under the assumption that color- kinematic duality exists. The main technical 
tool is a reformulation of the duality in terms of linear maps, modulo subtleties at loop level 
mostly inherent to the duality. This tool may have wider applications in both gauge and 
gravity theories, up to resolution of the subtleties. Here it is first applied to the large Britto- 
Cachazo-Feng-Witten (BCFW) shift behavior of gravity integrands constructed through the 
duality. Assuming color-kinematic duality and reasonable technical requirements hold these 
shifts are shown to be independent of loop order, which would imply massive cancellations 
with respect to the Feynman graph expression. More speculatively, the same approach is then 
applied to provide estimates of the overall degree of UV divergence in quite general gravity 
theories, assuming the duality exists. The cancellations obtained in these estimates depends 
on the exact implementation of the duality at loop level, especially on graph topology. Fi- 
nally, some evidence for the duality to all loop orders is provided from an analysis of BCFW 
shifts of gauge theory integrands through Feynman graphs. 
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1 Introduction 

Perhaps the largest unsolved problem in theoretical physics is the unification of the theory of 
general relativity with the principles of quantum mechanics. A range of possible solutions to 
this problem have been proposed. These typically involve either embedding into a much larger 
theory (e.g. string theory), appealing to non-perturbative dynamics (asymptotic safety) or 
by starting with a radically different formulation of the theory (loop quantum gravity), see 
e.g. [1] for a review. None of these have to date led to a phenomenologically viable and/or 
theoretically well-understood theory of quantum gravity. 
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It is well known that classically Einstein's gravity theory can be formulated as a La- 
grangian field theory, just as the standard model of particle physics can be. However, the 
latter can be turned into a well-defined quantum field theory by perturbative quantization. 
Divergent integrals which arise in this process in the standard model can be given a phys- 
ical meaning through renormalization: all divergences can systematically be absorbed into 
redefinitions of coupling constants and field normalizations. For renormalization to work the 
divergences must be under control and well understood. In particular for predictive power 
renormalization should involve only a finite number of coupling constants: this implies a finite 
number of measurements can fix these constants and all further measurements are predictions. 

The usual argument why gravity cannot be quantized according to the rules of quantum 
field theory involves a naive counting of powers of loop momenta. Around a flat background 
metric a perturbation theory in the coupling constant can be formulated just as it can in 
a gauge theory leading to Feynman graphs. When calculating say an n-point scattering 
amplitude at I loops one encounters complicated divergent integrals, roughly of the type 



where i is the number of massless propagators in this particular loop. The maximal number 
of loop momenta in the numerator, corresponding to the most divergent integrals in gravity 
theory, is k = 2i. In contrast, in Yang-Mills theories this number is k = i: the integrals 
appearing in gravity theories are therefore much more divergent in the ultraviolet "Zj —> oo" 
limit than their gauge theory equivalents. This difference can simply be traced to the fact 
that in four dimensional gauge theory the coupling constant is dimensionless, while it has a 
negative mass dimension in gravity in this number of dimensions. In particular this implies 
that in a typical calculation in gravity there are ever more divergent integrals appearing at 
each consecutive loop order, naively requiring an infinite amount of coupling constant and 
field renormalizations. If this number is truly infinite a perturbative quantum theory has 
no predictive power. In other words, there is no sensible implementation of renormalization 
known in gravity theories. There are two related ways out of this which preserve renormaliz- 
ability: the contributions of the divergent integrals could sum to zero or the divergences are 
not independent. In the latter scenario a symmetry would relate the divergent coefficients 1 . 

Investigating either scenario just mentioned is technically difficult in a Feynman graph 
approach in any gravity theory: there are many graphs, each with complicated expressions 
even for the simplest tree level calculations. At loop level this leads to complicated sums over 
divergent integrals. Using textbook methods it is therefore prohibitively difficult to find any 
would-be pattern within the divergences of gravity theories. Despite this, it has been verified 
non-zero divergences arise at one loop in generic matter-coupled gravity [2] and at two loops 
in pure Einstein gravity [3]. In supersymmetric theories of gravity no explicit divergence has 

An example of this are supersymmetric gauge theories on standard Af = 1 superspace: the field of mass 
dimension zero appearing in the vector multiplet can always be gauge transformed away provided the full 
supersymmetric gauge symmetry is present in the Lagrangian. 
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ever been demonstrated 2 . Beyond explicit finite order results powercounting in these theories 
remains therefore a rough upper bound on the apparent degree of ultraviolet divergence of 
gravitational theories: cancellations could very well be hidden within the sums. Until very 
recently however no mechanism was known to drive and study these would-be cancellations 
in loop diagrams. 

This has changed with the work of Bern, Carrasco and Johannson (BCJ) [11], [12]. 
These authors have proposed that at any loop order the integrand of gravity theories can 
in a precise sense be constructed as a double copy of two gauge theories. The field content 
of the gravity theory is simply the tensor product of the gauge theory factors. The driving 
force is a certain duality between color and kinematics which gives the kinematical factors the 
structure of a Lie algebra. The primary motivation for this type of construction at tree level 
are the Kawai-Lewellen-Tye (KLT) [13] relations in string theory where they are basically 
an observation about the relation between open and closed string correlation functions. The 
status of color- kinematic duality is reviewed in section 2. 

Since gauge theory amplitudes are well-behaved in the ultraviolet, a question arises if the 
BCJ representation of the gravity integrand as a double copy of gauge theory implies any 
generic cancellations with respect to naive powercounting. This question has been studied in 
specific example amplitudes, see footnote 2, up to relatively high (four) loop order. Although 
this progress is impressive it is also clear that the loop-by-loop approach of these examples 
will never lead to a proof of finiteness of any theory as this requires all order results. In 
this article a systematic study is initiated of the cancellations color-kinematic duality implies 
within quite general gravity theories to all loop orders in the gravity integrand. To this 
end we reformulate color-kinematic duality into a problem involving linear maps and study 
limits of this problem. The assumption that color-kinematic duality holds translates into the 
assumption a solution to this problem exists. The linear algebra approach initiated here could 
be of wider interest: in this paper it is used to study two examples of powercounting to all 
loop orders. The first is the study of large so-called Britto-Cachazo-Feng-Witten (BCFW) 
shifts of the gravity integrand, while the second makes some inroads into powercounting the 
overall degree of ultraviolet divergence directly. 

This article is structured as follows. First some of the necessary background material will 
be reviewed in section 2. In section 3 it will be shown that the improved scaling of tree level 
gravity amplitudes follows directly from color-kinematic duality. In section 4 it will then be 
argued that if color-kinematic duality holds at any loop order in addition to some technical 
assumptions, then the gravity integrand scales the same as the tree level amplitude. More 
speculatively, the same techniques will be applied to study powercounting aimed directly at 
UV divergences in section 5. In section 6 it will be shown how this scaling behavior can to 
some extent also be understood from Feynman graphs directly. The discussion and conclusion 

2 Recent results of explicit calculations include four point amplitudes in AT = 8 up to four loops [4], five 
point amplitudes up to two loops [5]. Beyond this there are firm arguments the first divergence appears at 
seven loops, see [6] and also [7]. In M = 4 explicit calculation has reached three loops, four points. See [8], [9], 
see also [10] for a more general argument. 



- 3 - 



section 7 follows. Some technical details are contained in the appendices. 



2 Review of concepts 

This section contains a brief overview over several concepts needed in the main body of this 
work. 

2.1 Color-kinematic duality 

Tree level closed string amplitudes in a flat background can be related to a certain sum over 
products of open string amplitudes. This is the content of the KLT relations [13], which basi- 
cally follow directly from the CFT picture of the string theory and the factorization of closed 
string vertex operators into left and right sectors. In the field theory limit these relations 
relate gravity amplitudes to a "square" of gauge theory amplitudes: a sum over products of 
gauge theory amplitudes with momentum-dependent constants. From a Lagrangian point of 
view these relations have always been quite mysterious. Recently, the relations have been 
recast in a more streamlined form by Bern, Carrasco and Johansson (BCJ) [11] at tree level 
and conjecturally extended to the loop level integrand [12]. In this form the relations are a 
consequence of a certain duality between color and kinematics. 

To expose color-kinematic duality amplitudes in a general gauge theory coupled to adjoint 
matter in D dimensions are rewritten as an expression in terms of cubic graphs only. Let 
Ti denote the set of all possible connected cubic graphs which can be drawn for a scattering 
amplitude with n external legs. Then for each graph in this set one constructs an associated 
color factor Cj by combining the structure constants of the three vertices in the obvious way. 
In this way the n-point Yang-Mills tree level amplitude can be written as, 



Here the coefficients rtj are referred to as kinematic numerators. Sj is the product of all 
propagators naturally associated to each cubic graph. This amounts to a rewriting of gauge 
theory amplitudes, as for instance in any gauge one could reabsorb quartic vertices into 
effective cubic ones according to their structure constant structure. 

The first non-trivial step in the Bern-Carrasco-Johansson approach is to demand that 
whenever color factors obey a Jacobi relation, the corresponding kinematic numerators must 
do so as well, i.e. 



It is useful to note that this identity can be given a nice graphical form depicted in figure 2.1. 
To state this first part of color-kinematic duality explicitly: 

Theorem 1 For every tree level amplitude in a theory with adjoint matter only one can find 
numerators rij such that they satisfy the Jacobi relations of equation (2.2) and reproduce the 
amplitudes through (2.1). 





Ck - Cj =^ Hi = Hk- Uj 
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Figure 1. A graphical illustration of the Jacobi relations. 



These numerators will be termed "color-dual" . The representation of a gauge theory ampli- 
tude using color-dual numerators in equation (2.1) will be called BCJ representation in the 
following. 

The numerators in equation (2.1) are not unique: for every set of shifts Aj such that 

Y^— = (2.3) 
r 

the numerators given by 

n'i = ni + Ai (2.4) 

describe the same amplitude through equation (2.1). If the shifts Aj also obey equation (2.2) 
the new numerators satisfy the Jacobi relations if the old ones did. This freedom in shifting 
numerators will be referred to as generalized gauge transformations. To some extent this 
shifting freedom is a manifestation of the fact that the color factors Cj are reducible: they 
obey (many) Jacobi relations. It will be shown below there is a bit more to it than this. The 
usual gauge freedom of amplitudes is a subset of the generalized gauge transformation. 

A generic set of numerators n\ which describe the amplitude through equation (2.1) will 
not be color-dual, i.e. 

Ci = Ck - Cj n\ - n' k + n'j = u m (2.5) 

will hold. However, in theories for which there are numerators nj which obey (2.2) there will 
be a generalized gauge transformation Aj for which 

Ac- 

Cj = c k - Cj Aj - A fc + Aj = u m ^2 = ( 2 - 6 ) 

holds. In this article when not explicitly stated otherwise generalized gauge transformations 
will be color-dual. 

The second step in color-kinematic duality is the following 

Theorem 2 Given a set of kinematic numerators which satisfy the kinematic Jacobi relations 
one can construct a gravity amplitude M by replacing the color factor in equation (2.1) by 
another set of kinematic numerators 

An = 9 n y - m 2 E—^ M n = - n ~ 2 E — ( 2 - 7 ) 
r s i r s i 

*■ i *■ i 
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with k the gravitational coupling constant. The field content of the gravitational theory is the 
tensor product of the field content of the two gauge theory copies. 

The generalized gauge invariance guarantees that shifts Aj of the kinematic numerators 
which satisfy (2.3) leave the squaring relations in equation (2.7) invariant. The reason is that 
the only input in (2.3) for the color factors is the algebraic Jacobi identity and the numerators 
are assumed to satisfy this identity. Note that the gauge invariance allows one to have only 
one set of numerators in equation (2.7) be explicitly color-dual. The other can be related to 
an explicitly color-dual set by a generalized gauge transformation as in equations (2.5) and 



The field content of the two copies of gauge theory does not have to be the same. For 
instance with both gauge theory numerators taken from J\f = 4 the BCJ construction gives 
J\f = 8 supergravity, while with numerators in pure Yang-Mills the resulting theory is that 
of a graviton, dilaton and a two-form. The latter theory is sometimes referred to as N = 
supergravity. Taking one set of numerators from pure Yang-Mills and the other from N = 4 
SYM gives M = 4 SUGRA. A more complete classification for D = 4 can be found in [14]. 

The two steps above are theorems since explicit sets of color-dual numerators have been 
found for instance in [15-17]. Furthermore, under the assumption that there are no non- 
trivial poles hiding in the set the double copy of equation (2.7) has been proven in [18]. A 
construction for numerators involving self-dual Yang-Mills theory has appeared in [19] as well 
as a construction in terms of auxiliary scalar field theories in [20] which will be discussed below. 
It has been shown at tree level for self-dual Yang-Mills theory that the symmetry structure 
implied by color-kinematic duality can be understood in terms of a hidden infinite kinematic 
Lie algebra [19, 20] which basically amounts to diffeomorphism invariance in a restricted part 
of space-time. Recently, there has also been work on a version of color-kinematic duality in 
three dimensions [21, 22]. Moreover, there has also been work on the color- kinematic duality 
involving polygons of MHV amplitudes [23]. 

More conjecturally the notion of color-kinematic duality has been extended to loops at 
the level of the integrand [12]. The starting point is as at tree level to express the integrand 
of a gauge theory amplitude at a fixed loop order I as a sum over trivalent graphs, 



Here D is the space-time dimension and Si indicates the symmetry factor of the trivalent 
graph i. The part under the integral sign will be referred to as the integrand of the gauge 
theory amplitudes. 

Note the color factors still obey all possible color-Jacobi relations. Then there is a natural 
extension of color-kinematic duality to the integrand in two steps. First 

Suspicion 1 For every loop level amplitude in a theory with adjoint matter only one can find 
numerators rij such that they satisfy the Jacobi relations of equation (2.2) and reproduce the 
integrand of the scattering amplitudes through (2.8). 



(2.6). 




(2.8) 
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Second, 

Suspicion 2 Given a set of kinematic numerators which satisfy the kinematic Jacobi rela- 
tions one can construct the integrand of a gravity amplitude M by replacing the color factor 
in equation (2.8) by another set of kinematic numerators 

M l n = ^ L f\{d D L^~ (2-9) 
J j=i r\ 

Assuming numerators rij can be found such that color kinematic duality (2.2) is satisfied 
unitarity implies that the double copy construction works as well, up to a subtlety identified 
below concerning generalized gauge transformations. This outline of a proof is again under 
the assumption that a set of color-dual numerators can be found which in addition has no 
spurious singularities. In fact, the existence of such a set can be taken as an additional, 
stronger, conjecture. Both assumptions are at loop level still an open question. Several finite, 
low multiplicity examples for amplitudes have been reviewed in footnote 2. New results on 
color-kinematic duality for low-point form factors up to four loops have recently been obtained 
in [24]. See also [25] for the integrand of the finite one-loop amplitudes in pure Yang-Mills. 

Note that it is natural to enforce a similar choice of loop momentum routing for the 
integrands on both sides of (2.8). This convention will usually be understood throughout 
this article. At loop level the numerators are not unique either. To quantify this freedom, 
suppose as before they are shifted by some amount Aj. In order for this shift to leave the 
amplitude invariant they have to fulfill the gauge condition 

i=i r, 

The shifts can be made to satisfy the Jacobi relations as well. The gauge condition is not as 
specific as above since the integrand only has to satisfy 

El A^c^ , x 

— = integrates to zero (2-11) 

1 i 

This leads to the question which terms vanish after integration. There are several ways terms 
can vanish after integration. The first is simple: the result might be a total derivative. An 
example of this is two terms that sum to zero after a shift in integration variable for one of 
them. I.e. . 

dL L2 (L + pi+p2 )2 " (L _ pi) 2 (L + p2 )2 =0 ( 2 - 12 ) 

Second, an integrand for a scattering amplitude can have no non-vanishing unitarity cuts. 
This is a standard reasoning to show this integrand integrates to zero. An example of this is 
the constant 1, which famously integrates to zero in dimensional regularization, 



d D L 1 = (2.13) 
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Vanishing of the cuts can be either purely algebraically or still involve a shift of the integration 
variable respecting the cut conditions. 

These varying ways of being able to satisfy the gauge condition (2.10) through vanishing 
terms have different implications for the squaring relations at the integrand level (2.10). 
Where at tree level it is clear that all color-dual shifts Aj that obey equation (2.3) will leave 
the squaring relation invariant, at loop level there is a problem with terms vanishing after 
integration since the numerators are generically momentum dependent. Hence the product 
of the gauge transformation with a kinematic numerators does not necessarily vanish. The 
example to keep in mind is the following type of integral, 

/oo rco 
xdx = while / x 2 dx ^ (2.14) 
-oo J — oo 

The safe terms are those which vanish algebraically on cuts which fix the loop momenta 
completely. In this context it is useful to note that it was found in [4] in order to make 
color-kinematics work in explicit loop computations, terms have to be added to the Yang- 
Mills integrand whose color-factor vanishes and whose integrand integrates to zero as well. 
These terms turn out to be crucial in order to obtain the correct gravity amplitude using 
the double-copy construction. If they were excluded the gravity integrand would not have all 
correct unitarity cuts. A full discussion of this freedom is necessary and interesting but is also 
outside the scope of this article. If one adds the further restriction that local numerators can 
always be found and these should be put into the double copy formula, then the ambiguity 
in the numerators is most probably to a large extent fixed, as it is in the known loop level 
examples. 

Recently, all- loop evidence for the conjecture was found investigating the IR structure 
of Yang-Mills and gravity [26]. Some more evidence from a direct analysis of BCFW shifts 
for the integrand is contained in section 6. In the following we will assume that color-dual 
kinematic numerators can be found at any loop order unless otherwise stated. Moreover, it 
will be assumed that for our results the subtlety in the role generalized gauge transformations 
in the double copy relations at loop level can safely be ignored. In particular, we assume that 
the double copy formula of equation (2.9) can simply be applied in this form to the numerators 
which will be obtained below. 

2.2 Solving Jacobi relations: the D 3 M basis 

To quantify the freedom in obtaining numerators in an equation like (2.1) it is useful to solve 
all Jacobi relations between the color factors. This amounts to considering the total system 
of Jacobi relations as a massive system of linear equations of the type 

Ac = AijG {-1,0,1} V i,j (2.15) 

with c a (2n — 5)!! dimensional vector of color structures. The kernel of the matrix A is a 
vector space which has several different baseses. One of these is particularly well-known [27] 
for the tree and one loop level and will be reviewed in this subsection. Obtaining an explicit 
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solution to the system at higher loop orders is generically an open problem, but for most 
purposes of this article only the existence of a solution will be needed. 

The system can be solved as follows: first single out two particles to be special, say 1 and 
n. Then a solution to the Jacobi relations follows by using them to maximize the number of 
structure constants encountered on the unique line through each trivalent graph from 1 to n. 
For this maximization one simply draws the trivalent graphs as the line connecting particles 1 
and n with all possible tree-graph subgraphs branching off from this line. These adornments 
can be 'dissolved' into the connecting line using the Jacobi relation in graphical form of figure 
2.1. The result is a single line with all external particles attached directly (see figure 2). 




Figure 2. Schematic application of the Jacobi relations (red ellipsis) to maximize the distance between 
two singled out legs (denoted by the hat) in a given graph. 

The resulting set is minimal since further application of any Jacobi identity would lead 
to non-trivial tree branches growing out of the connecting line. Hence the Jacobi relations 
can be solved in terms of the following set of structure constants, 

{/« 1 , Ma) « 1 / ai MS )« a • • • r-Vcn-W Vct e p ^ ■ ■ ■ ' n - x > ( 2 - 16 ) 

where P denotes the group of permutations of particles {2, . . . , n}. For future purposes define 

2, 3, . . . , U) = faia^aif 1 a 3 a 2 ■ ■ • f " 1 a n -ia n (2-17) 

This set has (n— 2)! different elements: all permutations of particles 2 through n—1. Scattering 
amplitudes can be expressed in this basis 

An = g 11 - 2 J2 F 0-M2), (7(3), . . . , cr(n - 1), n)A c °(l, a(2), a(3), . . . , a(n - 1), n) 

<76P{2,...,n-l} 

(2.18) 

where the coefficients A co were identified as the more well-known color-ordered amplitudes 
as was shown by Del Duca, Dixon, and Maltoni (D 3 M) in [27, 28]. This is straightforward 
to prove in a modern way using on-shell recursion. As shown in [27], the consistency of this 
equation with other possible choices of pairs of singled-out particles is a result of the Kleiss- 
Kuijf [29] relations for tree level color-ordered amplitudes. In fact, these relations can be 
seen as a consequence of the fact that a trace-based decomposition overcounts the number 
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of independent group theory structures available at tree level. Note that a similar argument 
holds at one loop where a minimal basis obtained this way is a ring of structure constants 
[27], while the corresponding relations for color-ordered amplitudes were found much earlier 
in [30]. The D 3 M basis at one loop for gluons, or more generally speaking particles in the 
adjoint, is explicitly given by 

A l = 1 =g n £ F(a(l),a(2),a(3),...,a(n)) / F°(a(l), a(n)) (2.19) 

aeP{l,...,n}/(Z n xZ 2 ) ^ ' 

where F is the ring of structure constants (sometimes also called adjoint trace) given by 

F(<r(l), a(2), a(3), . . . , a{n)) = n^lj^^-.J^H, (2-20) 

and I co denoting the single trace color-ordered integrand. The sum above runs over the 
permutations of the n external legs with inversions and reflections modded out. It is clear that 
in principle a solution to the Jacobi relations exists. An explicit solution would be interesting 
to find as this connects directly to extensions of the Kleiss-Kuijf relations to higher loop 
orders, see [31-33] for work in this direction. An alternative solution to the Jacobi identities 
minimizing the distance between two given particles is presented in appendix A. 

2.3 BCFW shifts and on-shell recursion 

In recent years on-shell recursion relations [34] [35] have been developed for gauge and gravity 
amplitudes at tree level and further steps to recursing loop level integrands and integrals 
have been taken, see [36] for a review. The main power behind on-shell recursion is that it 
reconstructs an amplitude from its residues at (a subset of) its kinematic poles. To do so a 
complex parameter z is introduced into the amplitude by shifting the momenta of any two 
legs of the amplitude, while keeping momentum conservation satisfied, by 

Pi=Pi + zq Pj=Pj~zq (2.21) 

The shift vector q is chosen such that it keeps the masses of the legs invariant, i.e. 

Pi-q = Pj-q = q-q = (2.22) 

giving two complex solutions for q. The amplitude is now a function of a complex parameter 
z and the original amplitude ^4(0) can be recovered by a contour integral around the origin 

Under the assumption that z = is an isolated singularity the contour of integration can be 
extended to infinity to give (neglecting the prefactor) 

A(0) = j> dz^- = - Yj (finite ^- (z = oo) (2.24) 

residues residue 
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The residues at tree level for finite z are products of lower point amplitudes summed over 
all internal states by unitarity. The residue at infinity lacks a similar physical interpretation, 
but if it can be shown to vanish (2.24) constitutes an on-shell recursion relation: the right 
hand side only involves tree amplitudes with a lesser number of legs. Whether a theory obeys 
recursion relations consequently depends on the behavior of the amplitude for z — > oo. If an 
amplitude falls off fast enough, i.e at least as ~ { z ~ l ), the residues at infinity are certainly 
absent and the amplitude can be computed through on-shell recursion. There has been work 
on establishing on-shell recursion even in the presence of boundary terms [37] and more 
recently there has been progress on on-shell recursion of Berends-Giele currents [38]. 

The large-z behavior of a Yang-Mills amplitude can be derived in several ways. All of 
them are variants of powercounting, i.e. tracing explicit powers of z in Feynman diagrams. 
In doing so the shortest (simply connected) path between two shifted legs along which the 
z dependence flows will be called "hard line". It can be shown that for a renormalizable 
gauge theory in D > 4 dimensions [39] that under a shift of two color- adjacent momenta a 
color-ordered amplitude scales as 



where /, are polynomial functions in {z~ l ) and fz^v antisymmetric in its indices. The vectors 
e are the z-dependent polarization vectors of the shifted legs i and i + 

The BCFW shift and on-shell recursion can be extended to the loop level at the level of 
the integrand [40], [41]. For the shift the above result in gauge theory follows from evaluating 
the amplitude through its Feynman graphs in the natural q^A^ = space cone gauge. The 
result for the shift for loop integrands is unchanged from the tree level result, i.e. it is given 
by (2.25). For integrands the notion of hard line follows from a particular choice for the 
routing of loop momenta. 

BCFW scaling analysis of gravity amplitudes at tree level and beyond 

The emergence of on-shell recursion relations for gauge theory amplitudes raised the natural 
question if similar constructions also hold for gravity theories. Here the outlook was originally 
very negative: simple powercounting of the involved Feynman graphs for BCFW shifts gives 
for n gravitons a z n ~ 2 divergence at z — > oo. Study of examples in [42] and [43] showed, 
however, that the recursion relations could actually hold, which was subsequently proven 
in [44] in four and [39] in D dimensions. See [45] for the extension to matter couplings. 
These works showed the actual divergence is a double copy of equation (2.25): the gravity 
polarization tensors can simply be decomposed in polarization tensors of two copies of gauge 
theory. A particular lucid argument for scaling of gravity amplitudes using a background field 
technique appeared in [39] where the scaling of an Einstein gravity amplitude was proven to 




(2.25) 
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be 



lim M(z)~ \e»(i)e»(i)] [e*(i + l)£*(f+l)l (z V ^f 1 {l/z) + z f 2 ^(l/z) + O(l/zf 

2— s>oo L JL J V / 

'zTHvMl/z) + z f 2t[ u>(l/z) + 0(1/*)) (2.26) 



where I? is as before an antisymmetric tensor. The product of polarization tensors in brackets 
belongs to the graviton, dilaton and two form. In passing we note that a form of the gravity 
tree level amplitudes which demonstrates fairly explicitly this scaling behavior has recently 
been found in four dimensions [46], [47]. 

This much-improved behavior under BCFW shifts at tree level compared to powercount- 
ing Feynman graphs immediately raises the question by what physical symmetry this is driven. 
Furthermore, based on experience with gauge theory integrands it is natural to investigate if 
the results obtained for tree level amplitudes can be carried over to the integrand level. It will 
be shown below, up to some assumptions, that the driving force of the cancellations in gravity 
is color-kinematic duality. Moreover, it will be argued the shift of the gravity integrand is 
indeed the same as the tree amplitude implying extensive cancellations taking place within 
the perturbative sum. 

2.4 Generalized inverses 

The concept of generalized inverses will be central in our study of BCFW shifts of color-dual 
kinematic numerators below. This is a well-developed mathematical subject of which some 
elements will be briefly reviewed here. See e.g. [48] for a more complete introduction. 

First generalized inverses for matrices are discussed. The purpose of introducing matrices 
is to solve systems of linear equations, 

Ax = b (2.27) 

Here A generically maps a vector space V to itself. If A is a square matrix of rank n it is well 
known that a unique inverse A -1 exists defined by 

A' 1 A = I = A A' 1 (2.28) 

which can be used to solve the equation to give 

x = A~ l b (2.29) 

In the case of singular or non-square matrices an inverse in this sense does not exist. 
However, the system of equations given by Ax = b can still have solutions: it is generically 
under- or overdetermined. This realization then leads to the theory of generalized inverses 
of matrices as was introduced first by Moore and later independently Penrose [49, 50]. Their 
work generalizes the notion of a matrix inverse to the cases of singular and rectangular 
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matrices, hence called 'generalized inverse'. The generalized inverse A + of a n x m matrix A 
is a m x n matrix defined here by the following condition only, 

AA + A = A (2.30) 

which reduces to the A + = A -1 if A is square and of full rank. In the more general case 
where A is singular or rectangular there is a solution to the system of equations 

Ax = b (2.31) 

only if consistency conditions hold. For instance, multiplying left hand side with a transpose 
vector c T for which c T A = holds should give c T b = 0. If A is a matrix these conditions can 
be shown to be equivalent to 

AA + b = b (2.32) 

This equation is certainly necessary as can be checked by multiplying left and right hand side 
in equation (2.31) by AA + . If the consistency conditions hold then the most general solution 
to the linear system is given by 

x = A + b + (I-A + A)y (2.33) 

with y an arbitrary n dimensional vector and I is the appropriate n x n identity. It is easy 
to check that (I — A + A)y is in the kernel of A. Less obvious is that it spans the kernel, 

kerA = {(I-A+A)y\y£C n } (2.34) 

Moreover, it is straightforward to check that (2.33) is a solution to the linear system since 

Ax = A{A + b + (I - A + A)y) = AA + b = b (2.35) 

The generalized inverse A + specified by equation (2.30) is not unique as it can be changed by 
transformations which leave equation (2.30) invariant. More precisely, if A + satisfies (2.30) 
then 

(A + )' = A + + (J - A + A) Y + W {I - A A + ) (2.36) 

also satisfies equation (2.30) for arbitrary matrices Y and W of size m x m and n x n 
respectively. Note that the extra pieces are related to the kernel of F. 

Example of a generalized inverse 

As an instructive example consider a singular matrix A given by 




(2.37) 



A generalized inverse A + of this matrix is given by 




(2.38) 
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as it obviously satisfies (2.30). The matrix problem 



Ax = b (2.39) 

can only have a solution in this case if 

AA + b = b -> b= (^\ (2.40) 
for some constant j3. The matrix I — A + A is simply 

I-A+A=(^fj (2.41) 

and its easy to see that 

keiA = {(I-A + A)y\yeC n } (2.42) 
indeed spans the kernel of A. Hence the most general solution to Ax = b of the form 

x = A + b + (/ - A + A)w (2.43) 

reads in this example 

where (5 was defined in equation (2.40) and w is arbitrary. The reader is encouraged to 
experiment with some alternative generalized inverses in this example generated by equation 
(2.36). 

The Moore-Penrose pseudo-inverse 

In this paper only the property (2.30) of the generalized inverse as introduced above is needed, 
as well as its ability to solve linear equations. In the literature more specific notions of 
generalized inverse abound which satisfy various other properties in addition to the defining 
equation (2.30). A prominent example of this is the Moore-Penrose pseudoinverse [49, 50]. 
This is a generalized inverse of a matrix A £ C mxri denoted by X which satisfies 

1. AX A = A (The definition of the generalized inverse) 

2. XAX = X 

3. (XA)* = AX (XA is hermitian) 

4. (AX)* = XA (AX is hermitian) 
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It can be shown [49, 50] that the solution to these conditions is unique. Other notions 
of generalized inverse impose a subset of these conditions. In practice the Moore-Penrose 
pseudoinverse can be computed using the singular value decomposition, i.e. one can factorize 
the matrix A as 

A = UT,V* (2.45) 

where U is a m x m complex unitary matrix, S is an m x n rectangular diagonal matrix with 
nonnegative real numbers on the diagonal, and V* is the conjugate transpose of a matrix V, 
i.e. an n x n real or complex unitary matrix. Then the pseudoinverse is given by 

X = VY, + U* (2.46) 

with E _1 the pseudoinverse of a, formed by replacing every nonzero diagonal entry by its 
reciprocal and then taking its transpose. Conveniently, Mathematica has the built-in functions 
SingularValueDecomposition and Pseudoinverse to calculate Moore-Penrose inverses. In the 
following only the notion of a generalized inverse as defined by (2.30) will be needed. 

Extension to more general linear maps 

In the above exposition a main assumption is that the matrix A maps a vector space V unto 
itself. The concept of generalized inverse however can be generalized beyond this to maps 
between two in principle very different vector spaces V\ and V%. In mathematics, this area is 
much less developed than the matrix case. 

To get an idea of what is involved, consider the space of integers V\ = Z and the space 
of integers modulo 2, V2 = Z mod 2. Let A be a linear map, 

A : Vi -»• V 2 (2.47) 

given by multiplication by 3, 

A(x) = 3x (2.48) 

Now consider the equation 

A{x) ~ b (2.49) 

Here the similarity sign is used to stress that this is an equation in the space V2. By construc- 
tion, b can only take the values or 1. In this case, it is easy to construct solutions directly 
for b either or 1: 

b ~ : x = 2k and if b ~ 1 : x = 2k + 1 for k G Z (2.50) 
It is not hard to recognize the kernel of the linear map A, 

kerA = 2k for k G Z (2.51) 
The generalized inverse A + in this case is a map V2 — > V± such that 

AA + A(x) ~ A(x) Vx G Z (2.52) 
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In this case, any 1-to-l map in the space V2 can be used as the generalized inverse. In this 
example the map A is invertible: there is a solution for every b. However, for the setup under 
study this is not always true. Consider 

C : V x -»• V 2 I C(x) = 2x (2.53) 

It is obvious that an equation of the type 

C{x) ~ b (2.54) 

for this map only has a solution if b ~ 0. For the generalized inverse again any 1-to-l map in 
the space V2 will do to satisfy equation (2.52). Note that also a map which maps any number 
unto will work. A consistency condition of the equation obtained by acting on left and right 
hand sides with CC + is 

CC + {b) ~ b (2.55) 

which is indeed only satisfied if b ~ 0. In this case this equation is apart from necessary also 
sufficient. The most general solution for the equation C{x) ~ b is 

x = C + b + ker C (2.56) 

In the example the kernel of the map A arises as the two spaces V\ and V2 have a different 
intrinsic dimensionality: the linear maps involved are in general many to one. This changes 
the setup slightly compared to the matrix case. Let A be a general map between vector spaces 
V\ and V2. Then the generalized inverse of A is the linear map A + for which 

AA + A{x) ~ A{x) Vx G Vx (2.57) 

for an equivalence relation on the space V2. The equation 

A{x) = b (2.58) 

has in general a solution only if 

AA + {b) ~ b (2.59) 

In general this is not sufficient as the maps involved are in general many to one. The more 
general consistency condition is that for every linear map C : V2 — > V2 such that 

C(Ax) ~ =>- C{b) ~ (2.60) 

A solution of equation (2.58) is given in terms of the generalized inverse by 

x = A + b + kerA (2.61) 

As before, the generalized inverse is not unique and can be modified by kernel mappings. 
More explicitly, if A + satisfies equation (2.57), then so does 

A + -> A + + Ci + C 2 (2.62) 
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where Ct are maps V2 —> V\ such that 



AC lX ~ Vx e y 2 (2.63) 

and 

C 2 Ax ~ Vie^i (2.64) 
hold. This is the linear map analog of equation (2.36). 

3 BCFW shifts of gravity amplitudes from gauge theory 

In this section the following question will be studied: can color-kinematic duality be used 
to derive BCFW shifts of gravity tree level amplitudes from gauge theory amplitudes? The 
main tool is the realization that the complete first step in the color-kinematic duality can 
be rephrased as a problem in linear algebra. This will lead to a singular system of linear 
equations which can be inverted in the generalized sense as introduced above to give the 
kinematic numerators n in terms of color-ordered amplitudes. By BCFW-shifting the large- z 
scaling of the numerators can then be read off and by using the double copy construction the 
large- z behavior of the gravity tree amplitudes will be obtained. In this section the kinematic 
numerators are always assumed to be color-dual. 

3.1 Warmup: four-point amplitudes revisited 

Let us illustrate the strategy at four points. The Yang-Mills four-point amplitude in the cubic 
BCJ representation is given by 

n s c s ntCt n u c u 
A A = + — + (3.1 
s t u 

where s, t, and u are the Mandelstam variables at four points and q the color factors given 
by 

The rii denote the corresponding kinematic numerators. Both, color factors and kinematic 
numerators satisfy the same Jacobi relation, i.e. 

c u = c s — c t and n u = n s - n t (3.3) 

Alternatively, the four-point amplitude can be written in terms of the D 3 M basis (2.18) 
singling out legs 1 and 4 

A 4 = c s A{\2U) + c u A(1324) (3.4) 

where A denotes color-ordered amplitudes. As these two representations are equivalent one 
obviously has 

nsCs + ntCt + nuCu _ + CuA ( 1324 ) (3.5) 

s t u 
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Inserting the Jacobi relations for color factors and kinematic numerators one the left hand 
side can be used to eliminate fit and q and arrives at 

c jus + ns^nu\ +C J^_ ns-nuX = CsA{12M) + CuA{im m 
\ s t J \ u t J 

which can nicely be written as a matrix equation 

jFij>n? = &Al => F^vP = Ai (3.7) 

with 

'-P4YA) -ft) 

i 7 is symmetric in its indices. Its determinant vanishes on-shell by momentum conservation 

det(F) = s + t + u = o (3.9) 
stu 

and has rank unity. For this reason the previously introduced concept of a generalized inverse 
(compare section 2.4) has to be applied to invert (3.7). One particularly simple representation 
is given by 

f+ = (o 4^) (3 - io) 

For a solution to the linear problem in equation (3.7) to exist the color ordered amplitudes 
have to satisfy 

FF + A = A (3.11) 
Spelling this out (using momentum conservation s + t + u = 0) one finds 



A(A(1234)-A(1324)»)+A(1324)^ ^ /a(1234)\ 
\^A(1324) + i('A(1234)-A(1324)f) J ~ \A(1324)J 



(3.12) 



This boils down to 

A(1234) - A(1324)- = (3.13) 
s 

which is nothing but the BCJ relation at four points. In other words, in this four point example 
the BCJ relations are a necessary and sufficient condition that the kinematic numerators 
satisfying the Jacobi relation exist. The most general solution to equation (3.7) is given by 

(:)- (s 4«) O) + (-■!)(:) 

for some vector w. The second term involving the vector w spans the kernel of F. Note 
that this kernel has a physical interpretation: it is the space of (color-dual) generalized gauge 
transformations. This can be demonstrated by evaluating equation (2.3) using the solution 
to the Jacobi identities employed here. 
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The behavior of the color-dual kinematic numerators under BCFW shifts can now be 
investigated in the four particle case. It is natural to shift particles 1 and 4, i.e. those kept 
fixed in the D 3 M basis. Under this shift the matrix entries of F in (3.8) scale as z° + 0(l/z). 
The generalized inverse F + scales as 

lim F + ~ ( ° ° , x V z° (° ° | + 0(l/z) (3.15) 

\0t + O(l/z)J \0 1J W ' y ' 

while 

lim (1 - F+F) ~ f j ^ W+Ofl/*) (3.16) 

Finally, the large- z scaling of the numerators is given by from (2.33) as 

lim m = lim (f + (z) A(z) + ker(F)(z)) ~ lim + f^M") + Kernel (3.17) 

Note that the kernel vanishes in (3.7) and does not contribute to the gravity amplitude in the 
double copy formula. That means that up to gauge transformations the numerators of the 
four point YM amplitude scale like color-ordered amplitudes adjacently shifted. The scaling 
of these objects is known from equation (2.25) and so the scaling of the kinematic numerators 
is given by 

lim ni ~ e(i)^(l)Jz V ^fi(l/z) + z #r) + 0(1/ z) (3.18) 

where e are the gluon polarization vectors, fi(l/z) is a function in \/z and B^ u is an antisym- 
metric matrix. By the double copy formula in equation (2.7) this behavior of the numerators 
up to generalized gauge transformations can be squared to give the shift of the corresponding 
gravity amplitude in M = SUGRA, 

lim M 4 ~ e(i) M e(i) K e(4),e(4) A 

z— >oo 

(zrTfiVz) + z°B^ v + 0(l/zj\ (zrj KX f(l/z) + z°B kX + C(l/z)) (3.19) 

The reasoning of this subsection will be extended in this section to all multiplicity at tree 
level and in the next section to the integrand. Before doing so the construction of numerators 
through linear algebra will be highlighted as this might have wider applications. 

3.2 Rephrasing tree level color-kinematic duality as linear algebra 

The first step is to solve the system of Jacobi relations. That is, there exists a basis for the 
color factors c 3 which in the tree level case is known to have dimension (n— 2)!. Moreover, any 
color factor of a tree level trivalent connected graph can be expressed as a linear combination 
of these, i.e. a rectangular matrix Aij exists for which 

Ci = A u c> b j = l,...,(n-2)! 1= l,...,(2n-5)!! (3.20) 
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At tree level, the D 3 M basis was reviewed above as an example of such a basis. Since this is 
just solving Jacobi relations this solution can be used in equation (2.1) for both color as well 
as kinematic part. This leads to 



Si 



which in turn leads to 



A n = g n y ^Y^n ] F k3 c k (3.22) 



with the symmetric (n — 2)! x (n — 2)! matrix F 



F jk = j (3.23) 

Moreover, the amplitude itself can be expressed in the color basis, 

A n = ^A n>j c> (3.24) 

3 

Hence if color-dual numerators exist then 

F jk h k = Anj (3.25) 

must hold for a symmetric matrix F. Note this equation (but not the inverse discussed below) 
also appears in [51]. Hence if F were invertible it would be proven numerators always exist 
and that they are unique. 

However, as shown above in the four point example the matrix F is in general singular. 
Moreover, we strongly suspect but have not been able to prove in general that the kernel of F 
is exactly the set of relations generated by the fundamental BCJ relations. This is basically 
a conjecture formulated in different ways in [11], see also [51]. Its essence is that the BCJ 
relations are the necessary and sufficient conditions for a color-dual representation to exist 
for scattering amplitudes. One can reason as follows: since F is symmetric, every vector in 
the null-space generates a relation for the vector of color ordered amplitudes by contracting 
left and right had side of the above equation. Using the D 3 M basis the set of amplitudes 
on the right hand side is the set of color-ordered amplitudes with two particles adjacently 
ordered. The most general set of relations for these color ordered tree amplitudes known 
and proven through other means is the collection of BCJ relations. Vice versa, every BCJ 
relation is a null vector of F. This makes the dimension of the kernel of the matrix F at least 
(n — 2)! — (n — 3)! and the rank of F maximally (n — 3)! under the assumption color- kinematic 
duality holds. We have checked this explicitly up to six points. Note that the kernel of F is 
the set of (color-dual) generalized gauge transformations, as noted above in the four particle 
example. 
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The matrix problem in equation (3.25) can be solved in terms of a generalized inverse 
F + if and only if the consistency condition 

FF + A = A (3.26) 

holds. As argued above, these conditions should be equivalent to the BCJ relations. Moreover, 
these conditions are equivalent to the existence of a set of Jacobi-satisfying numerators at 
tree level, i.e. theorem 1. As there are examples of such sets, the conditions must hold. Hence 
the system (3.25) can be solved to give 

n = F + A n + (1 - F + F)v (3.27) 

for an arbitrary (n — 2)!-dimensional vector v. Plugging this solution into the double copy 
formula in theorem 2 yields the most general form for the gravity amplitude in this language 

M n = n n - 2 {n b ) T Fh b = K n ~\A n ) T {F + + (/ - F + F) Y + W (I - F F+)) A n (3.28) 

as the natural consequence of the double copy construction in equation (2.7). Since the 
amplitudes satisfy the consistency condition this can for practical purposes by shortened to 

M n = K n - 2 (A n ) T (F+) A n (3.29) 

Note that one particular example of a matrix F + can be read off from this equation by 
comparing to the KLT relation obtained in [52]. In general it would be interesting to further 
explore the relation of the linear algebra approach to the momentum kernel of [15] especially 
through its link to string theory. 

Note that this analysis shows that most features of color-kinematic duality are captured 
by the features of the matrix F and its generalized inverse. The properties of these matrices 
are largely independent of the intricacies of Yang-Mills or gravity theories: they arise in any 
theory which has trivalent vertices made out of two Jacobi-satisfying structure constants. 
Other examples of these have been studied before [19, 20]. A particular class of these are 
scalar field theories, constructed to consist of a single massless scalar and only a three vertex 
which consists of the product of two structure constants satisfying the Jacobi relations. In 
other words, its perturbative series is designed to yield equation (2.1) directly. Note that this 
class involves specifying two 'gauge groups' and is therefore infinite dimensional. This class 
of quantum field theories will be referred to as 'trivalent scalar theories'. 

3.3 BCFW shifts of gravity amplitudes constructed by double copy 

The general setup just given can now be used to generalize the results at four points to 
arbitrary points at tree level. For concreteness use the D 3 M basis fixing legs 1 and n, to set 
up the problem, 

E^= E c a A(l,a(2),...,a(n-l),n) = Y / &iA (3.30) 
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where Tj runs over all cubic graphs and S n ^2 runs over all permutations of (2, n — 1). As 
demonstrated above, this reduces to 

FijW = Ai (3.31) 

which can then be solved in terms of a generalized inverse F + . As argued above, this equation 
is consistent. 

To obtain the large- z scaling of the numerators a BCFW-shift of particles 1 and n will 
be applied. It is clear how F scales under BCFW shifts: as ~ z . However, this does not 
automatically imply in general F + scales as z°. The canonical counter-example is 





(\ o o\ 




(\ o o\ 
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oio 




^0 J 




\o o o y 



(3.32) 



It is clear that the generalized inverse in this example is discontinuous in the limit e — > 0. 
More generally, it can be shown [48] that if the kernel of a matrix A has a smooth limit, 
then the generalized inverse has a smooth limit as well. For the pseudo-inverse this can be 
seen clearly using the singular value decomposition in equation (2.45) which basically reduces 
a general matrix to a problem of the same type as the example just studied. More general 
inverses are a kernel transformation of this analysis. 

To obtain the scaling of the generalized inverse F + BCFW shifts it will be useful to express 
the Yang-Mills amplitude via a basis of non-trivial eigenvectors of F. The key to doing this 
is the observation mentioned above that amplitudes in the class of trivalent scalar theories 
satisfy the consistency conditions (3.26) and therefore span the non-trivial eigenvectors of F. 
As argued above, the dimension of this space is (n — 3)!. Therefore for sufficiently rich choices 
of scalar field theories any solution to equation (3.26) can be expanded in terms of them, up 
to generalized gauge transformations, 

(n-3)! 
K 

where n ym are the Yang-Mills numerators and the numerators on the right hand side must 
not be in the kernel of F, 

FijW'K / (3.34) 

and be linearly independent. Equation (3.33) is equivalent to an observation made in [20]. 
Now it is known how F acts on the combination in equation (3.33) through equation (3.25), 

(n-3)! 

A jXM = F 3i n ym = £ a R (9?) (3.35) 

K 

where Q K are the amplitudes in the choice of trivalent scalar theories labelled by K. The 
construction can also be set up in the other direction, starting with the amplitudes. The 
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point of this expansion here is that for the trivalent scalar theories, 

h R = F+e R + (1 - F + F)v R (3.36) 

holds. In the class of trivalent scalar theories, it can now be shown that F + scales as ~ z°, 
up to a transformation in the kernel of F, i.e. a generalized gauge transformation. That is, 
the large- z scaling of F and F + is given by 

lim F = z° + 0(1/ z) 

Z ?°° , n (3-37) 
lim F = z + 0(1/ z) up to terms from ker(F) 

To show this, consider the D 3 M basis at tree level which singles out the BCFW shifted 
particles l,n. All scalar field theory amplitude coefficients in this basis scale as z° for this 
shift. The numerators in the trivalent scalar theories manifestly scale as z , up to generalized 
gauge transformation. Hence F + must scale as z°. 

Prom equation (3.35) it therefore follows that for Yang-Mills theories the essential infor- 
mation about scaling behavior of the numerators is captured by the coefficients a. In the case 
at hand the D 3 M basis is employed that singles out particles 1 and n. Shifting these amounts 
to a color adjacent shift for each component of the left hand side of equation (3.35). Since 
the scalar field theories manifestly scale as z , this shows that the scaling of the coefficients a 
follows immediately from the BCFW shift of color- adjacent gluons shown in equation (2.25) 

lim a R ~ e(i)^e(h) u (z^f R (l/z) + z°B R u ) + 0(1/ z) (3.38) 

Plugging this into equation (3.33) yields the result that the large- z behavior of the kinematic 
color-dual tree numerators at arbitrary multiplicity is 

lim rii ~ lim e(i),e(h) u (zri^f(l/z) + z°B» u ) + 0(1/ z) + kernel (3.39) 

z— >oo z— >oo V / i 

up to generalized gauge transformations. In writing this result we have used the fact that 
all numerators can be expressed as linear combinations of the basis numerators. These linear 
combinations only involve numbers. Note that strictly speaking this is an upper bound on 
the scaling: it could scale better. The result of equation (3.39) is sufficient to obtain the 
BCFW shift of the gravity amplitude through equation (2.7). 

It follows immediately from the numerators or equivalently from equation (3.28) that the 
(extended) Einstein gravity tree level amplitude at arbitrary multiplicity scales as two copies 
of Yang-Mills as in equation (2.26). This result itself was obtained through the background 
field method in [39]. What is new here is the explicit proof that the mechanism behind 
the large suppression in z scaling compared to the naive powercounting result ~ z n ~ 2 is 
color-kinematic duality. 
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Inclusion of renormalizable matter 

So far only gluonic matter has been considered and the question can be asked what happens 
if one adds fermionic and / or scalar matter in the adjoint representation to the gluons. The 
answer is simple: the only point were the field content played a role in the above argument 
is the shift of the tree level Yang-Mills amplitude in the D 3 M basis. Now from the generic 
arguments above it follows that this basis also exists for quite generic fermionic or scalar 
matter transforming in the adjoint of the gauge group. This and the associated KK relations 
can also be fleshed out [53]. The BCFW shift of adjacent gluons on tree amplitudes is given 
by equation (2.25) for renormalisable couplings of scalar and fermionic adjoint matter [45]. 
Hence the entire argument just constructed goes through for this class of theories, assuming 
color-kinematic duality holds: the gravity amplitudes constructed by double copy scale simply 
as the scaling of the double copy of the gauge theory components as in (2.26). Actually, also 
shifts of other particles than gluons can be straightforwardly be 'squared' this way. 

Improved scaling for permutation sums 

One simple application of the results of this section to tree level amplitudes follows from 
considering BCFW shifts of two gluons color- adjacent either to a permutation or a cyclic sum 
over external legs. Inspired by earlier work in [54] it was further fleshed out in [55], [56] that 
these shifts show improved shift behavior. A shift of a pair of color- non- adjacent particles is a 
particular example of this class of shifts. The improved shift behavior has since been proven 
in [57] and [58] using the BCJ relations for scattering amplitudes. Since the BCJ relations are 
a consequence of color-kinematic duality, it is not surprising that the same conclusion follows 
directly from the results obtained above on shifts of numerators. 

The point is that the particles within the permutation or cyclic sum have to couple to 
the hard line directly: a permutation sum on a current with more than one on-shell particle 
vanish by the photon decoupling relation. If the amplitude is now expressed in terms of the 
color ordered version of the BCJ representation, it is clear that every additional particle on 
the hard line will lead to a \ suppression from the 'propagator' type terms since only trivalent 
graphs are involved. By the results just derived the numerators scale as the Yang-Mills tree 
amplitude up to terms which do not contribute to the Yang-Mills amplitude. Hence the 
result for suppression conjectured in [55] [56] follow immediately. In a real sense, this is the 
color-kinematic counterpart of the proof in [57] and [58] . 

Further comments 

It should be obvious that the results of supersymmetric generalization of BCFW shifts consid- 
ered in [59] (see also [60]) for gravity amplitudes also simply follow from the above argument. 
Finally, nothing in this section depends crucially on selecting the D 3 M basis: any basis will 
do. Actually, from a formal point of view it is somewhat more natural not to solve the Jacobi's 
and simply treat Jacobi's and equation (2.1) as a massive set of linear equations which may 
be inverted using generalized inverses. A similar comment applies to the basis of trivalent 
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scalar theories used to study the scaling of the generalized inverse. Although above a choice 
of (n — 3)! different numerators not in the kernel was made, for the argument to go through 
it is not needed to really specify this number: one could take an over-complete basis. The 
expansion coefficients a are then not unique, but this can with a bit of tedious argument 
through generalized inverses be shown not to influence the result. 

4 BCFW shifts of gravity integrands from gauge theory 

Given the above results on BCFW scaling of Einstein gravity tree amplitudes through color- 
kinematic duality the next step is to extend the analysis to the integrand level. It will be 
argued that essentially all the steps of the tree level derivation go through, up to several new 
subtleties that will be identified. 

4.1 BCFW shifts of kinematic numerators using generalized inverses 

The integrand of an Z-loop n-point Yang Mills amplitude can be written in a cubic repre- 
sentation, just as the tree level amplitude can through equation (2.8), reproduced here for 
convenience: 

J j=1 r . % I 

As before, the color factors q in this representation obey a set of Jacobi relations and following 
suspicion 1 it will be assumed that a set of kinematic numerators rij can be found that obeys 
corresponding Jacobi relations. 

As before, these relations may be solved in terms of h different basis color factors and 
numerators, denoted by c and h. At one loop for instance, the D 3 M basis of equation (2.19) 
is an explicit example of such a solution with dimension h = (n — l)\/2. Inserting the basis 
into the cubic representation it can be rewritten analogous to the tree level case as 

A ~ / f[ d % £ jr™ -> ^ I II dL 3 F km n™ (4.1) 
j=i r\ * 1 3=i 

where k and m run from 1 to h and F is a matrix whose entries are sums over products 
of scalar propagators. The matrix F is not unique but depends column by column on the 
definition of loop momentum. These definitions have to be taken into account carefully when 
summing up the scalar propagators. The left hand side can also be expressed in the chosen 
color basis 




It can be useful to think of the integrand as defined at least in principle through Feynman 
graphs. A new feature at loop level is that one cannot extract the equation 

l YM = F jm h m (does not hold) (4.3) 
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as a consequence of equation (2.8). The main difference to tree level is that while the numer- 
ators live in the space of vectors of functions of external and loop momenta, 

/ fl(Pi,Li)\ 

n~ f 2 (p t ,U)j (4.4) 

the integrands live in the space of vectors of functions of external and loop momenta, identified 
up to terms which integrate to zero. That is, 

' fi(pi,Li) \ 

I- hiViM) (4.5) 

where two functions / and g are equivalent, / ~ g if 

I [J / Ud D L jg (4.6) 

The space of integrands is therefore in a real sense smaller than the space of numerators. As 
a linear map the matrix F therefore can have a non-trivial kernel. 

Some of the ambiguity has to do with routing the loop momenta on the left and right 
hand side. For convenience it is natural to involve the same choice on the two sides as will 
be enforced below. For a BCFW shift this implies on both sides there a hard line is chosen 
as the minimal length path between the shifted legs. In general the equation to solve is 

IJ M ~ F jm h m (4.7) 

or equivalently in the space on which the numerators live 

lJ M + lf nish = F jm n m (4.8) 

For a solution to this equation to exist there could be consistency conditions if there are left 
null eigenvectors of Fn. Note that 'null' in this sentence is up to terms which vanish after 
integration. Actually, at one loop we strongly suspect the relations found in [55] are the full 
set of consistency conditions for local numerators. At higher loops similar relations may exist 
as argued in [55, 56] on the basis of the improved behavior under non-adjacent BCFW shifts. 

In the linear map approach the assumption that color-dual numerators exists is simply 
the assumption that all consistency conditions are satisfied. This is the analogue of equation 
(3.26) at tree level. If this holds then by equation (2.61) the numerators at loop level are 
given by 

n = F + {I YM ) +kerF (4.9) 

In the tree level case the terms in the kernel of F were generalized gauge transformations 
which can be shown not to affect the squaring relation. As discussed above in section 2 
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the same does not hold for integrand numerators. This is a problem inherent to the color 
kinematic duality. It will be assumed this issue will not affect the outcome of the double copy 
relation for the integrand in a meaningful way in the following analysis. 

A further assumption is that the non-trivial eigenvectors of F are spanned by trivalent 
scalar theories. This amounts to the assumption that the space of solutions to the consistency 
conditions at loop level for color- kinematic duality is spanned by this class of theories. At one 
loop it is certainly true that both trivalent scalar theories as well as Yang-Mills theories obey 
the relation found in [55]. Moreover, both show a factor of - improvement for non- adjacent 
over adjacent shifts. In equations this assumption states that the left-side kernel of the map 
F 

left kerF = {m j \m i F ij hj = 0} (4.10) 

is the same set for numerators from the trivalent scalar theories as well as Yang-Mills theory 3 . 

With this assumption the scaling of the generalized inverse F + up to gauge transforma- 
tions can be studied as before by expanding the Yang-Mills numerators as 

hi = ct^hf (4-11) 

where 

Fjinf + (4.12) 

holds and the numerators are taken to be linearly independent. Here unbarred indices run 
from 1 to the number of color-ordered integrands independent under (suitably generalized 
loop level) KK relations, h, and barred indices from 1 to the number of independent color- 
ordered integrands under the consistency conditions (the analogs of the BCJ relation at tree 
level) at the integrand level. This gives an expansion for the integrand of Yang-Mills theory 
as 

(/ + / vanish )i = ajG/ (4.13) 

where ®f is the integrand of the Jth trivalent scalar field theory. As before, the action of 
any generalized inverse F + on ®f can be inferred from the general solution 

n K = F + e R + kerF (4.14) 

In trivalent scalar theories the left hand side scales as z° up to a generalized gauge trans- 
formation. The right hand side can be more complicated. In a typical choice of color basis 
the shift of two legs the elements of the integrand vector are linear combinations of color 
ordered integrands. At one loop for instance in the D 3 M basis a shift of two legs will be 
either adjacent or non-adjacent. Hence for the class of trivalent scalar theories 

(4.15) 




3 Some evidence for this will be provided later by an estimate of the kernel dimensions in both cases from 
BCFW shifts. 
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up to a generalized gauge transformation at one loop. Note the improvement in scaling of the 
integrand for non-adjacent shifts leads to a less well behaved generalized inverse F + . Now 
for the Yang-Mills integrand at one loop one immediately obtains 

«S = F+ljV ~ e(i) li e(n) v (zrTfO-/z) + z°B^ + 0(1/ z) + kernel (4.16) 

from the results for adjacent and non-adjacent shifts for loop level integrands. These results 
are obtained up to contribution which vanish after integration. If these would scale worse than 
the behavior indicated in the equation, then these would satisfy the consistency conditions 
themselves. Hence these are generalized gauge transformations. 

At higher loops the same argument goes through under an additional assumption: if one 
of the coefficients in the integrand vector scales shows improved BCFW scaling behavior, then 
the corresponding gauge theory coefficients must do so as well. Generically these coefficients 
are expected to either involve adjacent or non- adjacent shifts for which the shift behavior 
is known. Also, from the results of [31-33] there is a (conjectured) basis for 4, 5, 6 point 
integrands which is simply a subset of the full set of color-ordered integrands. This is enough 
to imply the assumption. 

Up to the assumptions the scaling of the numerators up to generalized gauge transfor- 
mations is given by 

lim n k ~ e(l)„e(h)J z V ^ f(l/z) + z°B^) k + 0(l/z) + kernel (4.17) 

2— >0O \ / 

4.2 BCFW shifts of gravity integrands constructed by double copy 

The double copy relation at loop level can be rewritten in terms of the linear map approach 
from equation (2.9). This reads 

< ~ I II d % £ |— ~ / II d % "''-U" J (4-18) 
J 3=1 r, * Si J j=i 

where h and h are sets of kinematic numerators from two copies of gauge theory. Knowing the 
large- z behavior of the numerators at any loop order (4.17) and of F it follows immediately 
that under a shift of particles 1 and n the large- z scaling of the n-graviton gravity integrand 
X n is given by a double copy of (4.17) 

lim Z n ~ £ (l)^(n) ! ,£(i)^(n) i >fzV^/(l/2) + z{jf v B^ + B^)+ 

V (4.19) 
z°(B^B^) + -(B^B^j) + 0{l/z 2 ) 

up to the assumptions. In particular, it is equivalent to the tree level result (2.26). Let us 
stress again that this result is under the assumption the generalized gauge transformation 
encountered along the way do not interfere with the double copy formula. 

Note this result displays even more dramatic cancellations in the powers of z than at 
tree level: at I loop for n gravitons the large-z scaling should naively be given by ~ z n-2+2i 
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but as one can observe from the above formula the largest power in z is quadratic and 
moreover independent of the loop order. As in the tree level case the driving force behind 
these cancellations is the assumed existence of color-kinematic duality. 

Comments 

The large- z scaling result above (4.19) for the M = integrand remains valid if scalars or 
fermions in the adjoint are included. This follows from the argument for adding matter at 
tree level in section 3. The needed results for scaling of Yang-Mills integrands can be found 
in [41, 56]. Similarly, the comments in section 3 on supersymmetric shifts and some minor 
variations also apply to the arguments presented here. 

Prom the scaling of the inverse of F in the integrand case one can estimate the number 
of independent integrands under the analog of the BCJ relation at loop level. This follows as 
the dimension of the kernel of F does not change in the limit: this is basically equivalent to 
the statement that F + (z) ~ z° up to a generalized gauge transformation. In the large z limit 
in a scalar theory it is straightforward to estimate the number of independent integrands for 
n external legs under the BCJ relations is the equal to the number of independent integrands 
under the KK relation for n— 1 particles. This follows from the fact that the leading diagram 
in the large z limit is a Feynman graph coupling to a one-leg off-shell current which involves 
n—2 on-shell particles. For these the analog of the KK relation holds, while no additional BCJ 
relation are expected for this off-shell object. Note in particular that both in the trivalent 
scalar theories as well as in Yang-Mills theory these results are basically the same. See also 
section 6 for the gauge theory analysis. This provides a cross-check of the expectation that 
the consistency conditions for the scalar theories are the same as those for the Yang-Mills 
theory. 

List of assumptions 

To be clear, the result in this section is subject to a number of assumptions, listed here: 

• color-dual numerators exist (in particular any consistency conditions on the Yang-Mills 
integrands are fulfilled) 

• the generalized gauge transformations involved do not influence the double copy rela- 
tions 

• the trivalent scalar theories span the set of integrands which satisfy the consistency 
conditions. 

• in the chosen color-basis the coefficients of the integrand for trivalent scalar theories 
scale the same or worse (higher powers in z) compared to the corresponding coefficient 
for the Yang-Mills integrand 
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5 Estimates of UV divergences assuming color-kinematic duality 

The main motivation for the study of cancellations in gravity theories is the UV behavior of 
the theory. The techniques based on linear maps employed above for the study of BCFW 
shifts can also be applied to the study of U V divergences through powercounting. Here it will 
be argued that color-kinematic duality leads directly to improved powercounting behavior 
of gravity integrands, subject to roughly the same assumptions as listed at the end of the 
previous section. Here we will only study the overall degree of divergence of graphs, ignoring 
sub and overlapping divergences. In other words, loosely speaking we study the limit 

li — y oo Vi (5-1) 

It should be noted a natural expectation present in much of the literature on divergences is 
that if a certain divergence is absent, a form of the integrand exists which makes this absence 
manifest under this scaling. If powercounting on the other hand yields a possible divergence, 
this could still vanish after integration and summation of all the terms. 

Note that it is known that the study of large BCFW shifts is related to the ultraviolet 
degrees of divergence, even beyond the fact that they are both based on powercounting. The 
canonical example of this is the good behavior of gravity tree amplitudes under BCFW shifts 
which is directly responsible for the absence of scalar triangle integrals at one loop in J\f = 8 
supergravity [60] (see also [61]). Hence by the results obtained above for tree amplitudes the 
absence of triangles at one loop for any number of points can by the results of this article be 
understood in terms of symmetry: it is a consequence of color-kinematic duality. Through 
unitarity, this argument factors through to the integrand at any loop order: there can be no 
triangle graphs, or at least none that can be isolated by a unitarity cut. In particular there 
can be no triangle divergences in either maximally supersymmetric Yang-Mills or gravity 
theory. 

For the purposes of this section one can re-write equation (2.9) as 




with Ii one of the copies of the gauge theory integrand, while the other copy yields the 
numerators h. 

At one loop it is well known that maximally supersymmetric Yang-Mills theory diverges in 
D = 8, while pure Yang-Mills theory is logarithmically divergent in D = 4. More generically, 
the usual decomposition of integrands of gauge theories with massless matter in terms of 
massless box, triangle and bubble integrals shows clearly the UV divergences in these theories 
at the bubbles 4 . 

At higher (> 1) loops the critical dimension where maximally supersymmetric Yang- Mills 
diverges logarithmically is [63], [64] 

D c =j + 4 (5.3) 
4 Intcrestingly, this relates sums over bubble coefficients to tree level amplitudes, see [60] and especially [62] . 
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This has been explicitly shown to be saturated up to including five loops [65] for general 
color and to six loops in the planar sector [66]. In particular in four dimensions this theory 
is UV finite. This translates to a power count of the integrand of amplitudes in maximally 
supersymmetric Yang-Mills theory of 



at loop order I. The number of propagator powers is that of a four point amplitude with the 
maximal number of loop propagators. 

In pure Yang-Mills the situation is different. Here naive powercounting of Feynman 
graphs gives a logarithmic divergence in two dimensions. This divergence cancels in dimen- 
sional regularization. This can be seen from the background field method: the functional 
form of the two-point correlator of two (background) gauge fields which contains the two 
dimensional divergence has to be proportional to 



by the background field gauge invariance. The two-dimensional divergence would give rise to 
the first term, while the second cannot arise with a two-dimensional divergence. Hence by 
gauge symmetry the two dimensional divergences must cancel in the sum and consequently 
generic Yang-Mills theories are logarithmically divergent in four dimensions. 

The main observation in this section is that the analysis of large BCFW shifts using color- 
kinematic duality of the previous section can be adapted to the limit of large loop momenta. 
More precisely, the inner product of all loop momenta with external momenta is taken to zero. 
As this is a well-defined limit, it should have a reflection for the kinematic numerators, up to 
the subtlety of the generalized gauge transformations. This limit will allow us to provide an 
estimate the overall degree of divergence of the corresponding gravity theory. 

The analysis in this section is subject to a number of assumptions, listed here: 

• color-dual numerators exist (in particular any consistency conditions on the Yang-Mills 
integrands are fulfilled) 

• the generalized gauge transformations involved do not influence the double copy rela- 



• the trivalent scalar theories span the set of integrands which satisfy the consistency 



• if the integrands of the trivalent scalar theories to be considered show improvement over 
the leading behavior the same holds for the gauge theory integrands. 

The main difference to BCFW shifts is that where in that case also tree level propagators 
are involved, here the tree propagators do not contribute to UV divergences as they do not 
involve loop momenta. With the above assumption the limit behavior of the linear map F + 




(5.4) 




(5.5) 



tions 



conditions. 
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can be studied in the class of trivalent scalar theories. Since the numerators of these theories 
do not depend on loop momenta, this determines the scaling of F + up to generalized gauge 
transformation from 

n R = F+ (/trivalent,^) + kerF ^ ^ 

for the K-th trivalent scalar theory. A subtlety arises here as the scaling of F + is now 
determined from the scaling of the different members of the minimal-under-KK-relations 
integrand basis /. As argued before, one version of this basis is likely to be a subset of all 
color-ordered integrands. One would now expect the multi-trace terms at least to have better 
powercounting than the planar ones. In pure Yang-Mills this is related to renormalizability: 
the UV divergences must be such that they can be reabsorbed into a counterterm proportional 
to the Yang-Mills action. Hence, these divergences must be proportional to tree amplitudes 
and in particular be single trace. 

If all the integrands in the trivalent scalar theory basis the integrands scale as ~ (L 2 )~ s 
for some integer 6. Then F + ~ {L 2 Y '. However, if the integrand vector scales partly as, say, 
(L 2 )~ Sl and partly as (L 2 )~ S2 , then the generalized inverse scales F + ~ (L 2 ) Sl or F + ~ (L 2 ) 52 , 
depending on the column number of F + . In schematic equations, 

if I ~ ( J then F+ ~ ( (L 2 ) Sl (L 2 ) &2 ) (5.7) 

The scaling of F + is as ever up to generalized gauge transformations. This behavior then 
determines the powercounting behavior of the numerators of the gauge theory by expanding 
them in the trivalent scalar theory basis. 

The next step is to decide which trivalent graphs to include in equation (2.8) or put 
differently, which graphs come with non-zero numerators. This is less innocuous than it 
sounds. For maximally supersymmetric Yang-Mills theory for instance it is extremely natural 
to exclude all triangle graphs. Apart from yielding a restricted set of graphs this also simplifies 
the Jacobi identities. In particular, this will admit a solution of the numerator Jacobis which 
is smaller than the one obtained for color Jacobis. In other words, there is a rectangular 
matrix such that 

hj = N jM n M (5.8) 

Here j ranges over all elements of the set of solutions to the color Jacobi relations. Note that 
in known examples, see e.g. [12], all Jacobi numerators can be expressed in terms of just one 
or two 'master' graphs. This limits the range of the index M, #M, to be << It would 
be very interesting to understand this from a group theory point of view. Many formulas 
given above can be streamlined in terms of the n basis, but this will mostly not be essential 
for this section. 

One loop 

At one loop, excluding triangle graphs gives for maximally supersymmetric Yang- Mills theory 

F + ~ (L 2 ) 4 A/" = 4,1 loop (5.9) 
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from the trivalent scalar theories, as all the coefficients in the integrand vector in the D 3 M 
basis contain at least a box. This gives 

n ym ~ (L 2 ) N = 4,1 loop (5.10) 

A very similar analysis for less-supersymmetric Yang-Mills yields 

F + ~ (L 2 ) 2 M < 4, 1 loop (5.11) 

which implies 

n ym ~ (L 2 ) N < 4, 1 loop (5.12) 

Combined with the double copy formula (5.2) this immediately gives the result that the double 
copy formula will generically scale as the best-behaved gauge theory copy. In particular, it 
implies the no-triangle property of M = 8 supergravity as a consequence of that of N = 4 
super Yang-Mills through the duality. Moreover, it correctly gives a logarithmic divergence 
in N = supergravity in four dimensions as well as the correct critical dimension (D=8) 
for N = 4 supergravity at one loop. In other words, at this loop order all known critical 
dimensions in the literature are reproduced, under the assumptions listed above. Of course, 
the analysis here does not include a check if the critical dimension is saturated. 

These results for scaling behavior are, as always, up to generalized gauge transformations. 
Moreover, if there is a vanishing integrand involved in the connection to the numerators, 
potential worse scaling behavior of this integrand can as in the BCFW shift case be absorbed 
into a generalized gauge transformation. 

Higher loops 

At higher loops the graph topologies appearing in the duality have in general a large effect 
on the generalized inverse of the linear map F, even apart from the issue sketched in (5.7) of 
differently scaling trivalent scalar theory integrands. For this consider N = 4 at four points. 
If one restricts to trivalent graphs with the maximal number of loop propagators, 

F + ~ (L 2 f l+1 M = 4, max # loop propagators (5.13) 

and consequently 

n ym ~ (L 2 )' -2 , J\f = 4, max # loop propagators (5.14) 

from I = 2 loops onwards. After squaring, this leads to a three loop estimate of the critical 
dimension of N = 8 which is too divergent compared to the known answer. Luckily, the 
answer can be seen from the known color-dual form of the three loop integrand in [12]: there 
are graphs such as on the left hand side of figure 3 which have one loop propagator less than 
the maximal number. Including these in the trivalent graphs and assuming all integrands 
scale homogeneously in the large loop momentum limit gives 

F + ~ (L 2 ) 3 ' M = 4, 1 loop propagator less (5.15) 
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and consequently 

n ym ~ (L 2 )' -3 , TV = 4, 1 loop propagator less (5.16) 




Figure 3. Contact graphs appearing at three loops (left) and at four loops (right). 

Similarly, at four loop graphs of the type on the right hand side of figure 3 appear in 
the explicit result for the four point integrand [67] . This kills of an additional power of L 2 , 
leading to 

F + ~ (L 2 ) 3 ^ 1 TV = 4, 2 loop propagators less (5.17) 

and consequently 

n ym ~ {L 2 ) l ~ 4 , TV = 4, 2 loop propagators less (5.18) 

This gives an estimate of the critical dimension of D c = 5.5 of TV = 8 supergravity at four 
loops. Note that up to and including four loops it is highly plausible by the above argument 
that the numerators in the higher point integrands obey the same scaling relations and that, 
consequently, the critical dimension for these integrands is the same. 

At five loops a problem begins to set in: if TV = 8 has to have the same critical dimension 
as TV = 4 5 , then for four points a graph with 3 tree propagators is necessary. Restricting to 
graphs without internal triangles, the only possibility for this is a tadpole. We have explicitly 
verified that four possible five-loop tadpole graphs made only out of box and pentagon type 
internal graphs exist using DiaGen [68] . If these tadpole graphs appears in the sum of equation 
(2.8), then the critical dimension by the power count estimate described here is the same as 
TV = 4 super Yang-Mills at five loops. If these tadpole graphs do not appear, this leads 
directly to an expected seven loop divergence in TV = 8 supergravity in four dimensions. If 
the five-loop tadpole appears and no further improvements occur then the above reasoning 
leads to the possibility that TV = 8 supergravity diverges in four dimensions at 8 loops. 

There are several ways to improve overall UV degree of divergence behavior. One is to 
include the possible four loop tadpole graph which consists of boxes only, see figure 4. If this 
graph appears, then the above reasoning estimates that TV = 8 supergravity diverges in four 
dimensions at 10 loops. If one relaxes the constraint to graphs without internal triangles, 

5 It should be stressed here it could still be true that the powercounting gives too pessimistic estimates. 
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Figure 4. Possible four loop tadpole graph made of boxes. 



there are no problems in rinding trivalent graphs to produce more cancellations. For J\f = 8 
to have the same power count as N = 4 SYM to all loop orders, one should include trivalent 
graphs which have the same UV power count as the M = 4 answer. Of course, the numerator 
factors here should be such that no triangle subgraph can be isolated by unitarity cuts in the 
gauge theory integrand. 

Powercounting in half-maximal supergravity 

The above methods can also be applied to study gravity integrands in gravitational theories 
with less than N = 8 supersymmetries. One well-studied theory is half-maximal supergravity. 
It was recently shown to diverge in D = 8 at one loop and D = 6 at two loops. Moreover, it 
was shown to be finite in D = 4 at three loops [8-10]. The gravity integrand of N = 4 can be 
obtained via the double copy construction as the product of Af = and M = 4 Yang-Mills in 
two different ways, 

Mi ~ J ^ d D L^=W=» ~ J n d D L] I^F+ =Q I N=0 (5.19) 



or 



< ~ / II dPLjI^n"- ~ / J! d D L^=°NF+ =4 I M=i (5.20) 

where N is the matrix out of equation (5.8) which connect N = 4 and TV = kinematic 
numerators. Note that the first formula gives a critical dimension of D = 8 by the result 
listed above, while the second gives D c = 4. Moreover, since the N = 4 SYM numerators are 
known to scale as (L 2 ) at least to four loops, the latter formula estimates M = 4 supergravity 
to be logarithmically divergent up to four loops. However, more cancellations can be hidden 
in the sums, just as they are at one loop. 

From the first formula it follows that the scaling of the pure Yang-Mills numerators sets 
the critical dimension of N = 4 supergravity compared to M = 4 super Yang-Mills. Naive 
expectation would be that these numerators scale as ~ (L 2 ) if all integrands in the color 
basis scale the same which would give J\f = 4 supergravity the same UV divergence as J\f = 8, 
while it is known that the 2 and 3 loop critical dimensions are 6 and ~ 4 respectively. At 
two loops, this translates into one power of L 2 . This can be argued to arise as follows: two 
loops is the first instance multi-trace terms appear which cannot related to the single trace 
term. In Yang-Mills these are expected to scale better in the large loop momentum limit by 
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renormalizability of the theory in 4 dimensions. In the basis of trivalent scalar theories, the 
leading contribution to the triple trace terms contains three loop propagators more: there 
may not be internal color triangle or bubble terms as they are proportional to either the color 
5 ab or f abc . Hence inverting this and taking into account the better scaling of Yang-Mills 
easily gives additional powers of L 2 . In fact, from the M = 4 supergravity result it is easy to 
infer that the triple trace integrand in pure Yang-Mills at two loops should scale two powers 
of L 2 better than the planar counterpart. 

A sketch of renormalizability in Yang-Mills 

Gravitational theories with less than M = 4 supersymmetry are expected to diverge in four 
dimensions. Moreover, generic Yang-Mills theory is logarithmically divergent but renormaliz- 
able in four dimensions. In this short subsection we sketch an outline of how renormalizability 
could work in M = 0. 




Figure 5. The bubble graph divergence at one loop 

Since at one loop the numerators scale as (L 2 ) , the only UV divergent graph in the 
color-kinematic dual integrand in four dimensions is a bubble presented in figure 5. Based on 
dimensional analysis, the divergence can only be proportional to p 2 , where p is the momentum 
flowing into the bubble. This cancels off one of the two \ propagators. Hence the UV 
divergence at one loop must by the results obtained here be proportional to a sum over 
trivalent graphs which looks like the color-dual representation of the tree level amplitude. 
This is the hallmark of renormalizability. Note that one can also conclude that for the 
quadratic Casimir for the kinematic algebra scales as C2 ~ p 2 + 0{(L 2 )~ l ). 

The same reasoning immediately extents to concatenated bubble graphs at higher loop 
orders. There can be more graphs which contribute at higher loops, see figure (6). Moreover, 
it is known the j3 function of Yang-Mills theories contains non-planar corrections [69], see 
also appendix B in [24] for a discussion. These graphs should also appear, with the value of 
the corresponding higher Casimirs for the kinematical algebra fixed, again, by dimensional 
analysis. The sum over all divergences should, at least in color space, reduce to a divergence 
times the corresponding tree amplitude. It would be interesting to pursue this further. 

Note that squaring the one-loop Yang-Mills result gave a logarithmic divergence in J\f = 
supergravity in four dimensions at one loop, as expected. Interestingly, taking the simple 
square of the Yang-Mills numerator scaling and analyzing the resulting graphs as above gives 
instead of a tree level-type graph with a \ propagator a contact-type graph. This is very 
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Figure 6. divergent trivalent graphs at two loops. 

much like what one would expect from a local counterterm different than the Einstein-Hilbert 
action. It would be interesting to see what the expression for this graph gives. 

From the 2 loop ~ (L 2 ) behavior of the pure Yang-Mills numerators extracted from the 
M = 4 supergravity result one would naively expect a quartic divergence at two loops in 
N = supergravity, driven by the triple trace parts of the Yang-Mills theory. However, since 
these arise in the non-planar sector of the theory, there is an argument to be made that the 
actual divergence is milder. It would be interesting to study this further. 

6 Towards an off-shell understanding of numerator scaling 

In this section it will be argued that the large- z behavior of the kinematic numerators at 
tree and integrand level can actually be seen from Feynman diagrams directly using stan- 
dard power counting. The key point will be to use color Jacobi relations to relate different 
contributions so that cancellations are achieved in the sum over Feynman graphs and the 
large- z scaling of the kinematic numerators found in the previous sections will be - at least 
conceptually - reproduced. 

6.1 Comparison to BCFW shifts directly via Feynman graphs 

To obtain the large- z shift of kinematic numerators from Feynman diagrams directly Feynman 
rules in the Feynman- 't Hooft gauge will be used and the gluon propagator will be put in the 
lightcone AHK gauge q ■ A = (see figure 6.1). An advantage of using this gauge is that if 
all the unshifted legs are left off-shell the scaling obtained holds at the level of the integrand 
to all loop orders. For brevity only the leading large- z contribution will be considered in the 
following but the subleading ones can be treated along the same lines. Remember that in our 
conventions the z dependence flows along the shortest path between the two shifted legs. 

Example: four point Feynman graphs 

At four points tree level there are only three cubic diagrams to consider: s, t, and u-channel 
with the four- vertex absorbed into the cubic graphs according to the color factors. In the 
BCJ representation this is written as 

c s n s cmt c u n u 
A 4 = + — + 6.1 
s t u 

In the following two particle momenta will be BCFW shifted and the other two legs will 
remain off-shell. In this way the analysis extends to the integrand. More precisely it extends 
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Figure 7. Fcynman rules used for power counting. The vertices are in Feynman-'t Hooft gauge and 
the propagator in AHK gauge. 
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Figure 8. s, t, and u channel cubic graphs used in the four point example. Hats denote the BCFW- 
shifted legs. 



to the Feynman graphs of the integrand with maximally one hard propagator in between the 
shifted legs. For graphs with more hard propagators the power counting will become more 
complicated as the number of graphs increases but in spirit it is similar. Hence, on the level 
of the integrand the leading large- z behavior for this class of diagrams is in a BC J form given 
by 

_^ c s n s cmt c u n u 

X 4 = + — + 6.2 
s t u 

with 

c u = c s - c t => n u = n s - n t (6.3) 
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Writing down the Feynman graphs and shifting particles 1 and 2 one can extract the large- z 
behavior of the kinematic numerators straightforwardly and it follows for the leading power 
in z for the three numerators (with indices a and p meaning to be contracted into appropriate 
currents) 



While these numerators satisfy the Jacobi relation above at leading order in z (this is special 
for four points) they do not yet have the form one would expect from the discussion in the 
previous sections, i.e z times a metric between the shifted legs. To get to this form one has to 
make use of the color Jacobi relations. It can be used to shift all terms not proportional to the 
metric from n u to nt and n s in (6.2). In the former numerator these will cancel with the other 
terms not proportional to jf 1 " and in the latter numerator these terms will become subleading 
in the large- z limit. This corresponds to shifting terms from the four-vertex between different 
kinematic channels. After some algebraic gymnastics one arrives at 



which nicely mirrors (leading part) of the result (3.39) obtained using pseudoinverses. It can 
be shown in a similar way that the subleading pieces are antisymmetric. 

Higher points and inclusion of matter 

In principle this procedure can be done in a similar fashion at higher points including the 
subleading parts. As the number of Jacobi relations increases quite rapidly with the number 
of points this becomes, however, more and more intricate. We could reproduce the scaling 
(3.39) in this direct Feynman graph approach up to including six points and our suspicion is 
that this can be done at any number of points. Furthermore, similar steps can be repeated 
for scalars and fermions in the adjoint and one arrives at the same result as before. This was 
also checked up to including six points. 

6.2 BCFW shifts versus color-kinematic duality 

The problem with the direct approach using Feynman diagrams is that the numerators ob- 
tained in this way do not satisfy the Jacobi relations beyond four points directly, i.e. they 




(6.4) 




(6.5) 



obey 



{rii — rij + = u a , a = 1, #Jacobis} 



(6.6) 
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where the right-hand side is non- vanishing and different for each Jacobi relation. One can now 
make use of gauge transformations on the numerators to bring them into a BCJ satisfying 
form, implementing a numerator shift 

rii — > rii + Aj with — — - = Mi (6-7) 

% 

For instance at tree level this works at follows: Following [11] there are (2n — 5)!! kinematic 
numerators / color factors at n points tree level. These obey ("~ 3 )^ 2n ~ 5 ) ! Jacobi relations 
which involve quite a lot of redundancy so that the number of non-redundant Jacobi relations 
is given by (2n — 5)!! — (n — 2)!. In other words (n — 2)! kinematic numerators / color factors 
are independent. Having constructed the numerators using Feynman diagrams one would 
have consequently obtained (2n — 5)!! — (n — 2)! non-redundant Jacobi relations that are not 
satisfied. 

{m - nj +n k = u a , a = 1, (2n - 5)H - (n- 2)!} (6.8) 

By enforcing a generalized gauge transformation on all (2n — 5)!! kinematic numerators the 
right hand side of the non-redundant Jacobi relations can be brought to zero by requiring (in 
addition to the gauge condition above) 

{Ai - Aj + A k = u a ,a = 1, (2n - 5)!! - (n - 2)!} (6.9) 

so that the non-vanishing right side of the kinematic Jacobi relations gets cancelled. Based 
on pure counting the generalized gauge condition would give (n — 2)! conditions on the shifts 
so that one the number of total equations for the shifts A is 

(2ra-5)!!-(n-2)! + (n-2)! =(2n-5)!! (6.10) 

#Jacobi relations to bring to zero # gauge conditions 

i.e. (2n — 5)!! equations for (2n — 5)!! shifts which means that in principle there should 
be a solution to this system of equation. This of course is hard to handle already at five 
points as the expressions involved become quite unhandy. Moreover, this raises the question 
if applying the generalized gauge transformations that bring the numerators into a BCJ 
satisfying form respects the large- z scaling? This question cannot be satisfyingly answered 
for arbitrary multiplicity using this approach but we have checked that (at least) up to five 
points shifts can be found and are such that they do not spoil the large- z behavior but in 
light of the discussion in the previous sections it is quite likely that such shifts exist for higher 
multiplicity as well. Of course, the same arguments go through at the level of the integrand. 

7 Discussion and conclusion 

The interaction of color- kinematic duality and powercounting has been studied in quite general 
gauge and gravity theories. The general analysis has been applied in two related but logically 
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distinct contexts: that of the large BCFW shift and that of large loop momenta. The main 
message of this article is that modulo the assumptions detailed above color-kinematic dual- 
ity implies massive cancellations with respect to naive powercounting arguments for gravity 
Feynman graphs. While at tree level the BCFW shift of the gravity amplitude was known to 
show this pattern, we add here the understanding that these cancellations are a consequence 
of an underlying symmetry principle. At loop level also the scaling of the gravity integrand 
is a new and, from a naive powercounting point of view, unexpected result. The extension to 
estimates of UV divergences presented above opens a window to a far more general approach 
than hitherto employed. 

In general the problem of proving color-kinematic duality remains open at loop level. 
Proving this is crucial to strengthen the validity of the results of this article. The linear 
map approach explored above might give a natural starting point for this. This analysis of 
color-kinematic duality should also provide a clear answer what role the vanishing terms in 
the integrand play. Moreover, it should include a discussion of which trivalent graphs to 
include in the construction. Both of these play an important role in the analysis of the UV 
divergences presented above. This obviously is also a direction for future research to decide 
the question if M = 8 is a finite quantum field theory. Our results above indicate that for 
J\f = 8 to have the same powercounting behavior as M = 4 SYM one should include tadpole 
graphs at 5 loops in N = 4. Moreover, there is a clear indication that graph topology is an 
important issue directly related to this question. 

There are many interesting questions to guide further research which lead from this 
work. For instance, it is easy to suspect but hard to make precise that the fall-off for gravity 
integrands obtained above implies the existence of on-shell recursion relations for the gravity 
integrand. At the very least, it implies the gravity integrand can be reconstructed from its 
single cut singularities, but unfortunately these have never been identified in terms of lower 
loop integrands except in the special case of N = 4 super Yang-Mills theory [40]. As an 
extension of this, one can ask the question if the \ fall-off for gravity amplitudes imply the 
existence of bonus relations for gravity integrands, along the same lines as [70] at tree level. 
Next, one could wonder if counterterms in theories which do have UV divergences have a 
non-trivial interaction with color-kinematic duality as suggested by the results of [71]. 

More generally, the techniques explored here have the potential to sharpen the analysis of 
perturbative quantum gravity considerably. At the very least they clearly suggest that massive 
cancellations take place in the gravity integrand, driven by the conjectured duality between 
color and kinematics. As such they should clearly be pushed further toward establishing 
concrete foundations. 
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Appendix A Solving Jacobi relations: minimizing distances 

In section 2 it was shown how the color Jacobi relations could be solved while maximizing 
the distance between two fixed labels. A natural second solution follows by using the Jacobi 
relation to minimize the number of structure constants encountered from 1 to n. To see how 
this works consider for instance the structure constants in (2.16). Now the Jacobi relation can 
be used to relate all different orderings in this set to one another, at the cost of introducing 
structure constants with less distance between 1 and n. It will be advantageous to aim for 
the natural minimal lexicographic ordering although this is not completely possible beyond 
four points. As a first step consider particle 2 and an arbitrary but fixed ordering of the 
remaining n — 3 particles. The Jacobi relation can be used to shuffle particle 2 to the minimal 
lexicographically ordered position, i.e. 

i 

F(l, ai, a 2 , • . . , a;, 2, . . . , n) = F(l, 2, ai, a 2 , • • ■ , n) - ^J^J F(l, ai, . . . , e, . . . , n)F(e, aj, 2) 

3=1 e 

(A.l) 

at the cost of introducing lesser length structure constant contractions. Now the would-be 
basis consists of one structure constant contraction of maximal length n — 3 and structure 
constants of lesser length. The same procedure can now be repeated at these lower length 
orders as well. There are now two complications. First, as mentioned before it is not possible 
to express these lesser length structure constants in terms of lexicographically ordered sets 
only. To see this consider for instance five points: there are in total six independent color 
structures but only five minimally lexicographically ordered ones. Secondly, there are now 
external sub-graphs attaching to the main line. The number of these are determined by 
applying Jacobi-relations in these subgraphs as well: for our purposes only the distance 
between the particles 1 and n will be important. 

Explicitly, consider structure constants contractions with a distance k between the special 
particles 1 and n on an n point amplitude. First consider the case where there are k legs 
with one and 1 leg with n — k — 2 legs. For the latter there are (n — k — 3)! ways to attach 
structure constants modulo the Jacobi identity. Then one can consider the case with k — 2 
legs with one and two legs with n — k particles in total. All different possibilities count. 

As a check that this second basis is minimal one can calculate its dimension. Counting 
the number of basis elements of a distance of k between the special particles 1 and n for 
an n-point amplitudes in the second basis can be mapped to the problem of counting the 
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number of possibilities of having a permutation of (n — 2) elements with k different cycles. 
This number is the unsigned Stirling number of the first kind, s(n — 2, k). For these numbers 



n-2 



J2s(n-2,k) = (n-2)! 



(A.2) 



k=l 



holds. Hence the second bases has the same dimension as the first one. 

Using either basis one can express the amplitude as a sum over basis elements /*. Then 
the right and left hand equation (2.1) as a sum over these 



Since the basis elements are independent, this constitutes (n — 2)! equations in (2n — 5)!! 
unknown numerators [29] which is always solvable. This shows explicitly there is a large 
degree of freedom in choosing the numerators. This is the freedom which is probed by the 
generalized gauge transformation. 

To put differently, the coefficients on the right hand side of equation (A. 3) are n — 2 
functions of the numerators. For any set of numerators for which (3i(n,j) = oti, equation (2.8) 
gives the tree level amplitude A. 
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